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Abstract 



We study the relativistic Vlasov-Maxwell system which describes large systems of particles 
interacting by means of their collectively generated forces. If the speed of light c is considered 
as a parameter then it is known that in the Newtonian limit c — > oo the Vlasov-Poisson system 
is obtained. In this paper we determine the next order approximate system, which in the case 
of individual particles usually is called the Darwin approximation. 
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Plt 1 Introduction and Main Results 

43 , 

The relativistic Vlasov-Maxwell system 



( d t f + v ■ VJ + (E + c~ l v x B) ■ V v f = 0, 

cV x E = -d t B, cVxB = d t E + Anj, f . 

V-E = 47rp, V-B = 0, (RVMc) 

P ■= Jfdv, j := Jvfdv, 

describes the time evolution of a single-species system of particles (with mass and charge normalized 
to unity) which interact by means of their collectively generated forces. The distribution of the 
large number of particles in configuration space is modelled through the non-negative density 
function f(x,v,t), depending on position x G M 3 , momentum v G M 3 , and time t G R, whereas 

v = (l + c- 2 v 2 y 1/2 v e R 3 (l.l) 

is the relativistic velocity associated to v. The Lorentz force E + c~ x v x B realizes the coupling 
of the Maxwell fields E(x,t) G R 3 and B(x,t) G R 3 to the Vlasov equation, and conversely the 
density function / enters the field equations via the scalar charge density p(x, t) and the current 
density j(x,t) G R 3 , which act as source terms for the Maxwell equations. It is supposed that 
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collisions in the system are sufficiently rare so that they can be neglected. The parameter c denotes 
the speed of light, and J always means L 3 . At time t — 0, the initial data 

f(x,v,0) = f°(x,v), E(x,0) = E°(x), and B(x, 0) = B°(x) 

are prescribed. In this work we treat the speed of light as a parameter and study the behavior 
of the system as c —>■ oo. Conditions will be establish under which the solutions of (jRVMc|) 
converge to a solution of an effective system. We recall that in [21] it has been shown that 
as c — > oo the solutions of ()RVMc|) approach a solution of the Vlasov-Poisson system with the 
rate 0(c _1 ); see jU E] for similar results and [E] for the case of two spatial dimensions. The 
respective Newtonian limits of other related systems are derived in [2m H|- It is the goal of 
this paper to replace the Vlasov-Poisson system by another effective equation to achieve higher 
order convergence and a more precise approximation. This will lead to an effective system whose 
solution stays as close as 0(c~ 3 ) to a solution of the full Vlasov-Maxwell system, if the initial data 
are matched appropriately. In the context of individual particles, this post-Newtonian order of 
approximation is usually called the Darwin order, see [23 EH] and the references therein. Let us 
also mention that weak convergence properties of other kinds of Darwin approximations for the 
Vlasov-Maxwell system were studied in [SJ [2]. In the present paper we mainly view the Darwin 
approximation as a rigorous intermediate step towards the next order, where in analogy to the 
case of individual particles ^3] radiation effects are expected to play a role for the first time. Since 
at the radiation order the corresponding dynamics of the Vlasov-Maxwell system most likely will 
have to be restricted to a center manifold-like domain in the infinite dimensional space of densities 
(to avoid "run- away" -type solutions [23 E]), it is clear that several new mathematical difficulties 
will have to be surrounded in this next step. Then the ultimate goal would be to determine the 
effective equation for the Vlasov-Maxwell system on the center manifold, which should finally lead 
to a slightly dissipative Vlasov-like equation, free of "run-away" solutions; see [TTJ H2] for a model 
of this equation and more motivation. 

Compared to systems of coupled individual particles, for the Vlasov-Maxwell system one im- 
mediately encounters the problem that so far in general only the existence of local solutions is 
known. These solutions are global under additional conditions, for instance if a suitable a priori 
bound on the velocities is available; see the pioneering work [Hj, and also PHIO], where this result 
is reproved by different methods. This means that from the onset we will have to restrict ourselves 
to solutions of (|RVMc|) which are defined on some time interval [0, T] that may be very small. On 
the other hand, in [21] it has been shown that such a time interval can be found which is uniform 
in c > 1, so it seems reasonable to accept this restriction. 

In order to find the desired higher-order effective system, we formally expand all quantities 
arising in (|RVMc|) in powers of c _1 : 



/ 


= fo + c- 


Vi + c" 2 / 2 + • • 


E 


= Eo + c 


~ X E X + C 2 E 2 + 


B 


= B Q + c 


- 1 B 1 + c- 2 B 2 + 


P 


= Po + C 


Vl + C~ 2 p 2 + • 


J 


= 3o + c~ 


X Ji + c~ 2 j 2 + . . 



where Pk = J fk dv and jk = J vfk dv for k — 0, 1, 2, — Moreover, v = v — (c~ 2 /2)v 2 v + ... by 
(jl.lj) . where v 2 = \v\ 2 . The expansions can be substituted into (jRVMcJ) . and comparing coefficients 
at every order gives a sequence of equations for these coefficients. 
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At zeroth order we obtain 

Vx£ = 0, V-E = 47rp , Vx5 =0, V • B = 0. (1.2) 

If we set 

B = 0, 

then the Vlasov-Poisson system 

dtfo + v ■ V x / + E ■ V v f = 0, 
Eq(x, t) = — J \z\~ 2 z p (x + z, t) dz, 

Po = / /o dv, 
/o(x,v,0) = / (x J t;), 



(VP) 



is found, with z = \z\ l z. 

Next we consider the equations at first order in c -1 . Here 

V x Ex = -d t B = 0, V-£7i=47rpi, V x Sj = + 4vrjo, V • B x = 0, (1.3) 

needs to be satisfied for the fields; also see [T2] . Using (jl.2j) . we get Ai?! = — 47rV x j and therefore 
define 

Bi(x,t)= / |x-y| _1 Vx j Q (y,t)dy= / |^|~ 2 ^ x j (x + z, t) d^. (1.4) 



Regarding the density fi, we obtain the linear Vlasov equation 

dtfi + v ■ V x fx + Ex ■ VJo + E ■ V„/i = 0. 
Hence if we suppose that fx(x, v, 0) = 0, then we can set 

fx = and E x = 

consistently. 

The field equations at the order c~ 2 are 

V x E 2 = -d t Bx, V • E 2 = 4vrp 2 , Vx5 2 = + 47rji =0, V ■ B 2 = 0. 

Therefore we can define 

B 2 = 0. 

Calculating the equation for the density / 2 and taking into account (jl.3j) . we arrive at the following 
inhomogeneous linearized Vlasov-Poisson system, for which we choose homogeneous initial data: 

dth + v ■ V x f 2 -\v 2 v VJo + E ■ V v f 2 + (E 2 + vx B x ) ■ V J = 0, 

AE 2 = d 2 E + 4vr(Vp 2 + d t j ), (LVP) 
f 2 (x,v,0) = 0. 

At this point we need to discuss the solvability of the Poisson equation for E 2 . Restricting our 
attention to initial data f° for (jRVMc)) with compact support, it will turn out below that both 
p 2 and jo have compact support and thus lead to unproblematic sources. The first term d 2 E has 
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to be examined more closely. Since po(-,t) has compact support for all t, see (|2.1|) below, we can 
calculate the iterated Poisson integrals 



A-\d?E )(x,t) = -l-J^-d 2 t E (y,t) 



1 f dw . ., 

— / | — r / \z\ z o t po{x + w + z,t) dz 



Air J \w 
1 



4tt 

If dy 
2 
1 
2 



dyd 2 p (y,t) j du\y — x — u\ 1 \u\ 3 u 
J \yl x \ ( y ~ x "> % ''Po^) = \ J zdfp (x + z,t)dz (1.5) 




\z\ 1 (v — (z ■ v)z) d t f (x + z, v, t) dzdv, 



where we used (TVPJ) . eq. (|5.27|l from the appendix, and d t po + V ■ jo = in conjunction with an 
integration by parts, the continuity equation itself being a direct consequence of ()VP|) . In view of 
fTTKjl and ijTWjl we thus define 

E 2 (x,t) = ~^J z dt po(x + z,t) dz — J \z\~ 1 d t jo{x + z, t) dz — J \z\~ 2 z p 2 (x + z,t) dz. (1.6) 

By (jl.6|) . ()VP|) . and a further integration by parts, we obtain the alternative expression 

i r f. , o 



E 2 (x } t) = — j j \z\ z (3(z ■ v) —v ) fo(x + z,v,t) dz dv 



\z\~\l + z ® z) (E p )(x + z,t)dz - / \z\~ 2 zp 2 (x + z,t) dz. (1.7) 



The first aim of this paper is to show that 



/" := fo + c- 2 f\ 



E D 
B D 



2) 

2; 



E + c~ 2 E 2 , (1.8) 



yields a higher-order pointwise approximation of (jRVMcjl than the Vlasov-Poisson system; we 
call (jl.8)l the Darwin approximation. It is clear that for achieving this improved approximation 
property also the initial data of (jRVMc)) have to be matched appropriately by the data for the 
Darwin system. For a prescribed initial density f°, we are able to calculate (fo,Eo), Bi, and 
(f 2 , E 2 ) according to what has been outlined above. We then consider (|RVMc|) with initial data 

f(x,v,0) = f°{x,v), 

E(x,0) = E°(x) := E Q (x,0) + c~ 2 E 2 (x,0), (IC) 
B(x,0) = B°(x) := c~ 1 B 1 ( y x,0). 

Before we formulate our main theorem let us recall that solutions of (jRVMcjl with initial data (jlCj) 
exist at least on some time interval [0, T] which is independent of c > 1; see |21| Thm. 1], and 
cf. Proposition 12 .21 below for a more precise statement. This time interval [0, T] is fixed throughout 
the paper. 
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Theorem 1.1 Assume that f° G C°°(M 3 x M 3 ) is nonnegative and has compact support. From f° 
calculate (f ,E ), B\, and ($2^2), and then define initial data for hRVMc}) by 1/6]) . Let (f,E,B) 
denote the solution of \RVMc}) with initial data HIQ) and let (f D , E D , B D ) be defined as in M.ty) . 
Then there exists a constant M > 0, and also for every R > there is Mr > 0, such that 

\f(x,v,t)-f D (x,v,t)\ < Mc~ 3 (xel 3 ), 

\E(x,t) -E D {x,t)\ < M RC - 3 (|ar|<i2), (1.9) 
\B{x,t) -B D (x,t)\ < Mc~ 3 (xGM 3 ), 

for allv G IR 3 ; t G [0,T] ; and c > 1. 

The constants M and Mr are independent of c > 1, but do depend on the initial data. Note 
that if (jRVMc|) is compared to the Vlasov-Poisson system (|VPJ) only, one obtains the estimate 
\f(x,v,t) - f (x,v,t)\ + \E(x,t) - E (x,t)\ + \B(x,t)\ < Mc" 1 ; see ^ Thm. 2B]. 

Approximate models have the big advantage that, since by now the Vlasov-Poisson system is 
well understood, the existence of (fo,E ), and here also of B\ and (f2,E 2 ), does no longer pose 
serious problems; note that in (jLVP|) the equation for f2 is linear. Therefore one can hope to get 
more information on (jRVMc)) by studying the approximate equations. As a drawback of the above 
hierarchy, one has to deal with two densities / 0) fi an d two electric fields E , E 2 to define f D and 
E D . Therefore it is natural to look for a model which can be written down using only one density 
and one field. It turns out that the appropriate (Hamiltonian) system is 

f d,f + (1 - ic-V)t> ■ V,/ + (E + c-h< x B) ■ V„/ = 0, 

C \l E= 7 d f'. V '£ = 4 *". (DVMc) 
cAB = — 47rV x j, v ' 

K p = Jfdv, j = f(l - \ c~ 2 v 2 )v f dv, 

which we call the Darwin- Vlasov- Maxwell system. We note that (f D , E D , B D ) solves (jDVMc|) up 
to an error of the order c -3 . 



Theorem 1.2 Assume that f° G C°°(1R 3 x M 3 ) is nonnegative and has compact support. Then 
there exist c* > 1 and T* > such that the following holds for c> c* . 

(a) If there is a local solution of WVMc\l . then the initial data E° and B° of IIP VMc}) at t = 
are uniquely determined by the initial density f°. 

(b) The system hDVMc]) has a unique C 2 -solution (f*,E*,B*) on [0,T*] attaining that initial 
data (f°,E°,B°) att — 0. This solution conserves the energy 

n= j J (\v 2 - l - c - 2 v A ^tdxdv + ^ J (|V0*| 2 + |VA A*\ 2 ^j dx, 

where the potentials <fi* and A* are chosen in such a way that B* = V A A*, V • A* = 0, and 
-V0* = E* + c~ l d t A* . 

(c) Let (f,E,B) denote the solution of \RVMc}) with initial data (f°,E°,B°). Then there exists 
a constant M > 0, and also for every R > there is Mr > 0, such that 

\f(x,v,t)-f*{x,v,t)\ < Mc~ 3 (x G M 3 ), 
\E{x,t) - E*{x,t)\ < Mrc' 3 (|ar|<i2), 
\B(x,t) - B*(x,t)\ < Mc' 3 (i£l 3 ), 

for allv G 1R 3 ; t G [0, min{T, T*}], and c > c*. 



5 



Instead of performing the limit c — > oo in (|RVMc|) it is possible to reformulate Theorem 11.11 
in terms of a suitable dimensionless parameter. Taking this viewpoint means that we consider 
flRVMc|) at a fixed c (say c = 1) by rescaling a prescribed nonnegative initial density f°, for which 
we suppose that f° e C 00 ^ 3 x M 3 ) has compact support. To be more precise, let 

vf°(x, v) dx dv, 

where v is taken for c = e _1//2 ; cf. (11.1)1 . Then v is viewed as an average velocity of the system. 
Now we introduce 

f £ '°(x,v) = e 3/2 f (ex,e~ 1/2 v) 
and consider f £, ° for c = 1. It follows that 




v s 





vf e,0 (x,v)dxdv = \/e / / wf°(y,w) dydw = y/ev, 



i.e., the system with initial distribution function f £, ° has small velocities compared to the system 
associated to f°. Starting from /°, we next determine (fo,E ), Bi, and (f 2 ,E 2 ), and then the 
initial data for (jRVMcj) via (jlCj) with c = e _1//2 , as in Theorem ll.il Next we note that (/, E, B) is 
a solution of (|RVMc|) with c = e~ 1//2 if and only if 

f £ (x,v,t) = e 3/2 f(e Xl e- 1/2 v } e 3/2 t) } 
E £ (x,t) = e 2 E(ex,e 3/2 t), 
B £ (x,t) = e 2 B(ex,e 3/2 t), 

is a solution of (jRVMcj) with c = 1. We further introduce 



]{x, 


v,t) 


= e^foiex^-^v, 


e 3/2 t), 




E%l 


[x,t) 


= e 2 E (ex,e 3/2 t), 






B{i 


>,t) 


= e^B^ex^'H), 






\(x. 


v,t) 


= e"l 2 f 2 {ex,e-y\ 






E £ 2 ( 


[x,t) 


= e 3 E 2 {ex,e 3/2 t), 






Po< 


[x,t) 


= J fo(%,v,t) dv = 


- e 3 p (ex 


,e 3/2 t), 


Jo • 


[x,t) 


= J vfo(x,v,t)dv 


= ^ 7/2 Jo( 


ex,e 3/2 t) 




[x,t) 


= J /lO^M) dv = 


4 / 

= e p 2 \ex 


,e 3 > 2 t). 



Straightforward calculations then confirm the following statements: 

(a) (/o, Eq) is a solution to (|VPJ) with initial data f° if and only if (/q , Eq) is a solution to (|VPJ) 
with initial data f £, °, 

(b) Bi solves ABi = —4ttV x j if and only if B\ solves ABf = — 47rV x j'q, 

(c) (/ 2 , E 2 ) is a solution to (jLVPj) if and only if (/f , J5§) is a solution to 
ft/I + u • V,/! - ^ 2 « • V x f £ Q + £ £ • V„/f + {El + vxBf). VJ^ = 0, 

AE £ = d 2 E £ + 4vr(Vp 2 + d t f ), 
fl(x,v,0) = 0. 
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Therefore Theorem 11.11 mav be reformulated in a way which parallels |13( Thm. 2.2], where the 
case of individual particles is considered which are far apart (of order 0(e~ 1 )) and have small 
velocities (of order 0(y/e)) initially. Note that in this result the Lorentz force is determined up to 
an error of order 0(e 7 ^ 2 ), and the dynamics of the full and the effective system can be compared 
over long times of order 0(e~ 3 ^ 2 ); see [EU p. 448]. 

Theorem 1.3 Assume that f° G C°°(IR 3 x M 3 ) is nonnegative and has compact support. From f° 
calculate (/o, Eq), B\, and (/2, E2), and then define initial data for \RVMc\) by 1/6]) with c = e^ 1 ^ 2 . 
Let (/, E, B) denote the solution of \RVMc]) on [0, T] fore = e^ 1 ^ 2 with initial data IIG]) . Moreover, 
let f £) ° , f £ , E £ , B e , /o, Eq, B{, /|, E\, p £ Q , Jq, and p\ be defined as above. Then (f £ ,E £ ,B £ ) is a 
solution of \RVMc\) on [0, e~ 3 ^ 2 T] for c = 1 with initial data 

(A E £ , B £ )(x, v, 0) = {f<°{x, v), E £ (x, 0) + E\{x, 0), Bf(x, 0)). 

In addition, there exists a constant M > 0, and also for every R > there is Mr > 0, such that 

\f £ (x,v,t)-f £ (x,v,t)-fl(x,v,t)\ < Me 3 (16I 3 ), 

\E £ (x,t) - E £ (x,t) - E £ 2 (x,t)\ < Mrs 1 t2 {\x\<e~ l R), 
\B e (x,t) -Bl(x,t)\ < Me 112 (i6l 3 ), 

for all v G M 3 , t G [0, e~ 3 ^ 2 T], and e < 1. The constants are independent of e. 

By definition of the rescaled fields, these fields are slowly varying in their space and time variables, 
which means that we are considering an adiabatic limit. It is clear that also Theorem 11.21 could be 
restated in an analogous e-dependent version. 

The paper is organized as follows. Some facts concerning (jVPj) . (jLVP)) . and (jRVMc)) are 
collected in Section 121 The proof of Theorem 11.11 is elaborated in Section El whereas Section 0] 
contains the proof of Theorem 11.21 For the proofs we will mostly rely on suitable representation 
formulas for the fields (refined versions of those used in [HI 121]), which are derived in the appendix, 
Section EJ 

Notation: B(0, R) denotes the closed ball in M. 3 with center at x = or v = and radius R > 0. 
The usual L°°-norm of a function (p = (p(x) over 1 6 R 3 is written as \\<f\\ x , and if if = tp(x,v), we 
modify this to IM^,,. For m G N the W m ' oc -norms are denoted by |M| ma ., etc. If T > is fixed, 
then we write 

g(x,v,t,c) = O cpt (c- m ), 
if for all R > there is a constant M = Mr > such that 

\g(x,v,t,c)\ < Mc~ m (1.10) 

for |x| < R, v G M 3 , t G [0, T], and c > 1. Similarly, we write 

g(x } v,t,c) = 0(c~ m ), 

if there is a constant M > such that (jl.K)j) holds for all x, v G M 3 , t G [0,T], and c > 1. In 
general, generic constants are denoted by M. 
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2 Some properties of (EE), (EVE), and ( jkVMcl 



There is a vast literature on (jVP|) . see e.g. Sect. 4] or JH] and the references therein. For our 
purposes we collect a few well known facts about classical solutions of (|VPj) . 

Proposition 2.1 Assume that f° G C°°(K 3 x R 3 ) is nonnegative and has compact support. Then 
there exists a unique global ^-solution (fo,Eo) of HVP]) . and there are nondecreasing continuous 
functions Pvp, Kvp '■ [0, oo[— > R stzc/i that 

ll/oWiu<nriu, 

supp/ ( v ,t) C 5(0,iVp(t)) x B(0,P VP (t)), (2.1) 

II/o(*)II w + II^)(*)IIm<^vp(0. 

/or t G [0, oo[. 

This result was first established by Pfaffelmoser [TS], and simplified versions of the proof were 
obtained by Schaeffer and Horst j^j; a proof along different lines is due to Lions and Perthame 

For our approximation scheme we also need bounds on higher derivatives of the solution. This 
point was elaborated in [TH], where it was shown that if f° G C fc (R 3 x R 3 ), then (f ,E ) posses 
continuous partial derivatives w.r.t. x and v up to order k. The existence of continuous time- 
derivatives then follows from the Vlasov equation. Thus (fo,Eo) are C°°, if f° is C°°, and by a 
redefinition of Kyp we can assume that 

\\fo(t)\\ 3tXiV <K VP (t), te[0,oo[. (2.2) 

The existence of a unique C 1 -solution (f2,E 2 ) of ()LVP|) follows by a contraction argument, 
but we omit the details. Furthermore it can be shown that there are nondecreasing continuous 
functions Plvp, Klvp '■ [0, oof— > R such that 

supp/ 2 (.,.,t) C B(0,P LVP (t)) x B(0,P LVP (t)), (2.3) 
Wf2(t)\\ 1>x>v + \\E 2 (t)\\^<K LVP (t), (2.4) 

for t G [0, oof. 

Concerning solutions of (jRVMc)) . we have from |2"T) Thm. 1] the following 

Proposition 2.2 Assume that f° G (^^(IR 3 x M 3 ) is nonnegative and has compact support. If 
E° and B° are defined by \1C)) . then there exits T > (independent of c) such that for all c > 1 
the system \RVMc]) with initial data HIC)) has a unique C x -solution (f,E,B) on the time interval 
[0,T]. In addition, there are nondecreasing continuous functions (independent of c) Pvm,Kvm '■ 
[0,T] -»• R such that 

f(x,v,t) = if \v\>P VM (t), (2.5) 
\E(x,t)\ + \B(x,t)\ <K VM {t), (2.6) 

for allx G R 3 ; t G [0,T], and c > 1. 

In fact E 10 and B° do not depend on c in |2*Ti Thm. 1], but an inspection of the proof shows 
that the assertions remain valid for initial fields defined by (jl(J|) . 
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3 Proof of Theorem 11.11 

In Section 15 . 1 . II below we will show that the approximate electric field E D from ()1.8|) admits the 
following representation: 

E D = Eg t + E£ + E° d + O cpt (c- 3 ), (3.1) 

with 

E^ t (x,t) = — \z\~ 2 z(po + c~ 2 p 2 )(x + z, t)dz 

J\z\>ct 

—c~ 2 I \z\~ l d t jv(x + z,t) dz + - cT 2 / zd 2 p (x + z,t)dz, (3.2) 

J\z\>ct ^ J\z\>ct 

EP t (x,t) = -j \z\~ 2 z(p + c~ 2 p 2 )(x + z,i(z))dz 

— 2{z ■ v)z) fo(x + z, v, t(z)) dv dz 
\z\~ 2 (2(z ■ v)v + v 2 z — 3z(z ■ v) 2 ) f (x + z, v, t(z)) dv dz 
\z\~ (z (g> z — l)E p (x + z, t(z)) dz, (3.3) 
E° d (x,t) = c-\ct)- 1 [ f(z-v)zf°(x + z,v)dvds{z) 

J\z\=ct J 

+c~ 2 (ct)~ 1 / / {{z ■ v)v - (z ■ vfz) f°(x + z, v) dv ds(z), 



\z\<ct 

f 

\z\<ct ■ 



+c- 2 



+c~ 2 



'\z\<ct 

I 

l\z\<ct 



'\z\=ct 



where the subscripts 'ext', 'int', and 'bd' refer to the exterior, interior, and boundary integration 
in z. We also recall that z = \z\~ 1 z and t(z) = t — c~ 1 \z\. On the other hand, according to Section 
l5.1.2l below we have 

E = E cxt + E int + E hd + 0(c~ 3 ), (3.4) 



with 



E ext (x,t) = - [ \z\ 2 z (p + tdtpo + 7;t 2 d 2 po) (x + z, 0) dz 

J\z\>ct V I J 

z d 2 po(x + z,0) dz — c~ 2 / \z\~ l d t j®{x + z, 0) dz, (3.5) 



'\z\>ct 



J\z\>ct J \z\>ct 

E- m t(x,t) = — \z\~ 2 zp(x + z, t(z)) dz 



+c 



:\<ct 
-1 



+c" 2 



'\z\<ct 



'\z\<ct 



\z\ 2 / (2{z ■ v)z — v) f(x + z,v,i(z)) dv dz 



+c' 2 \z\-~ 1 (z®z- \){Ef)(x + z,v,i(z))dv dz, (3.6) 

J\z\<ct J 



E hd (x,t) = E° d (x,t). 



In order to verify (jl.9|) . we start by comparing the exterior fields. Let a; G 5(0, R) with R > 
be fixed. Then we obtain from ()3.5|) and ()3.2|) . due to \z\ = 1, and taking into account 



p 2 (x, 0) = / f 2 (x, v, 0) dv = 

by (1LW1) . as well as (Q, (Q, and Q, 

|£ ext (x, t)-E^(x,t)\ 

/ |z|~ 2 p (z + z, t) - p (x + z, 0) - td t p (x + z, 0) - - t 2 d 2 p (x + z, 0) 

J|z|>ct ^ 

+c~ 2 / |z|~ 2 |p 2 (^ + - P2(^ + 2,0)1 <fc 

J\z\>ct 

+c~ 2 / / \v\\d t fa(x + z,v,Q) — d t fo(x + z,v,t)\dvdz 

J\z\>ct J 



< 



dz 



+\c 2 f [ \d 2 f {x + z,v,t) - d 2 f {x + z,v,0)\dvdz 

Z J\z\>ctJ 

< M [ \z\~ 2 ( [ (t-s) 2 P VP {s) 3 K V p(s)l B{0tPvp{s)) (x + z)ds)dz 

J\z\>ct \Jo J 

+Mc~ 2 [ \z\~ 2 ( [ P LV p(s) 3 K LV p(s)l B ^p LVp{s)) (x + z) d^\ dz 

J\z\>ct \Jo / 

+Mc- 2 [ \z\~ X ( ( Pvp(s) 4 K VP (s)l B{0t p vp{s)) (x + z)d S S ]dz 

J\z\>ct \J0 / 

+Mc~ 2 [ ( [ P VP (s) 3 K V p(s)l B{0t p Vp{s)) (x + z)ds)dz 

J\z\>ct \J0 J 

< Mt 3 \z\- 2 l mR+MQ )(z) dz + Mtc' 2 \z\-\\z\- 1 + l + \z\)l B (o,R+M )(z) dz 

J\z\>ct J\z\>ct 
-3. 



< M R c- 3 ; (3.7) 
note that here we have used 

M = max (Pyp(s) + K VP {s) + P LV p{s) + K LVP {s)^ < oo, (3.8) 

and for instance 

t 3 [ \z\- 2 l mR+Mo) {z) dz < {ct)-H 3 [ \z\ dz < M R c~ 3 . 

J\z\>ct J\z\<R+M 

To bound \Ei nt {x,t) — E® t (x,t)\, we first recall from Thm. 2B] that 

\E(x,t) - E (x,t)\ =0(c- 1 ). (3.9) 

Actually the initial conditions in |2I] are different, but we only added terms of order c -2 , so that 
an inspection of the proof in j2I] leads to (J3.9j) . Next we define 

H(t) = sup{\f(x,v,s) - f D {x,v,s)\ : x G R 3 , v G IR 3 , s G [0,t]}, 
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as well as 



Mi = max (Pvm(s) + Pvp(s) + Plvp(s)) < oo. 
se[o,T] V / 



Then f(x, v, s) = f {x, v, s) = f 2 (x, v, s) = for x G M 3 , |v| > M u and s G [0, T). Also if i? > 
is chosen such that f°(x,v) = for \x\ > Ro, introducing the constant 

M 2 = Ro + TM l + max (p vp (s) + P LV p(s)) < oo 
se[o,T] V / 

it follows that f(x,v,s) = f (x,v,s) = f 2 (x,v,s) = for |x| > M 2 , v G M 3 , and s G [0, T]. Let 
x G 5(0, i2) with R > be fixed. From ipTBjl . Q , (I3~SJ) . (Q, flEJ), and < t(z) < t for \z\ < ct 
we obtain 



\E int (x,t) - EP t (x,t)\ 



< 



+c 



z\<ct 

1 



(/ - fo - c 2 / 2 )(£ + ^, u, *(z)) dv 



dz 



z\<ct 



+c 



+c 



+c 



+c 



(2(z ■ v)z -v)(f-f - c- 2 f 2 ) (x + z, v, t(z)) dv 



(2(z ■ v)z — v) f 2 (x + z, v, t(z)) dv 



dz 



(v 2 z + 2(z ■ v)v - 3(z ■ v) 2 z) (f-fo- c~ 2 f 2 )(x + z, v, t(z)) dv 



z\<ct J 

\ wHy 

z\<ct J 

\ wHy 

z\<ct J 

' (1 - z ® z)([E - E ]f)(x + z, v, t(z)) dv 
< M{Ml + Mt)H{t) [ \z\~ 2 1 b{0) m 2 ){x + z) dz 

J\z\<ct 



dz 



dz 



(v z + 2(z ■ v)v — 3(z ■ v) z) f 2 (x + z, v, t(z)) dv 



dz 



dz 



• .U.U ,'.!/„ c~ 3 
+MMlH(t) c~ 2 

• MM:; M (l c -4 



\z\<ct 

+MM 3 1 \\n\ XtV c- 3 

< M R (c- 3 + H(t)), 
since for instance 



/ \z\ 2 l B(0jM2) (x + z)dz 

~ 2 1b(o,m 2 )(x + z) dz 
\ \z\~ 2 1 B ( ,M 2 )(X + z) dz 

J\z\<ct 



\z\<ct 



1-1 



\z\<ct 



lB(0,Ah)( X + Z ) dz 



(3.10) 



/ \z\ 2 1b(o,m 2 )(x + z) dz < \ 

J\z\<ct J\z 



\z\~ 2 dz < Mr. 



z\<R+M 2 



Recalling that the Ebi(x,t) = E® d (x,t), we can summarize (|3.4jl . (|3.1|) . (|3.7|) . and f)3. lOj) as 

|£(x, t) - E D (x, t) | < M R (c~ 3 + (t)), (3.11) 
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for \x\ < R and t G [0,T]. Formulas (jSHJ), (pTT3|) . (|5^ . (PTT7|) . (ETTHj) . and lETTO!) . and an 

analogous (actually more simple) calculation also leads to 



\B(x,t)-B u (x,t)\ <M( C - S + H(t)), (3.12) 

for x G R 3 and t G [0,T]. It remains to estimate h = f - f D . Using (IRVMcjl . (jVPj), and 

flLVPl. it is found that 

d t h + v ■ V x h + (E + c^v x B) ■ V„/i 
= -dj 13 - v • VJ C - (E + c~ l v x B) ■ V J D 

= (y - - c~V u - vj ■ VJ + c~ 2 (t> - v) ■ VJ 2 

+(B D - E) ■ V v f + c-\E D - E) ■ V v f 2 - c-*E 2 ■ V„/ 2 

+c- 2 ((t; - v) x B x ) ■ V„/ + c-\v x (B D - B)) ■ V„/o - c" 3 ({) x B) • V„/ 2 . 

If |w| < Mi, then also \v\ = (1 + c~ 2 v 2 )~ l l 2 \v \ < \v\ < Mi uniformly in c, and hence 

v-(l -^c~ 2 v 2 ^v < Mc~ 4 . 

Next we note the straightforward estimate \Bi(x,t)\ < M for \x\ < M 2 and t G [0, T], with Bi 
from (H3J). In view of the bounds (|2~TJl . (Q, and JZEJ, thus by (gHH) and (|3~T2) . 

|<9t/i(x, u, t) + • V x h(x, v, t) + (B(x, t) + c _1 i) x B(s, £)) ■ V v h(x, v, t)\ 

< M(c~ 3 + H(t)) (3.13) 

for |x| < M 2 , \v\ < Mi, and t G [0,T]. But in {(x,w,t) : \x\ > M 2 } U {(x,v,t) : |u| > Mi} we have 
h = f — f D = by the above definition of Mi > and M 2 > 0. Accordingly, ()3.13|) is satisfied for 
all x G M 3 , v G M 3 , and t G [0,T]. Since h(x,v, 0) = 0, the argument from p. 416] yields 



H(t) < [ M(c~ 3 + H(s))ds, 
Jo 



and therefore H{t) < Mc~ 3 for t G [0,7]. Then due to (EHTl) and (l3~T2l . |B(x,t) - B D (x,t)| < 
Mr c" 3 for < R and t G [0, 7), as well as |B(x, t) - B D (x, t)\ < Mc~ 3 for x G M 3 and t G [0, T]. 
This completes the proof of Theorem 11.11 □ 



4 Proof of Theorem 11.21 

In this section we will be sketchy and omit many details, since the proof is more or less a repetition 
of what has been said before. First let us assume that there is a C 2 -solution (/*, E*, B*) of (jDVMc|) . 
existing on a time interval [0, 7*] for some 7* > 0, such that supp /*(•, •, t) C R 3 x R 3 is compact 
for all t G [0, 7*\. Then 



B*(x,t) = cT 1 J \z\- 2 zxj*(x + z,t)dz, (4.1) 
AE*(x,t) = 47rVp*(x,t) + c- 1 d t V x B*(x,t). (4.2) 
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Since f*(x,v,0) = f°(x,v) and j*{x,0) = j{l-\c~ 2 v 2 )vf*(x,v,Qi)dv = c- 2 v 2 )vf°{x, v ) dv, 

it follows that B*(x,0) is determined by /°. In order to compute the Poisson integral for E*, we 
calculate by means of the transformation y = w — z, dy = dw, and using ()5.27|) below, 

c^A-^dtV x B*){x,t) 

= -s? / v » x ( / i*r 2 * x + *■ 4 > *) 

1 



47TC 2 

1 f dw 

2d 2 



dw I dz\z\- 2 \x-w + z\- 1 [zV-(d t j*)(w,t)-(z-V)d t f(w,t) 
x-w]V- (d t f)(w,t) -([x-w]- V)d t f{w,t) 



\x — w\ 



1 f dz . , „ „ . , 

= y z®z)d tJ *(x+z,t). 

If we invoke the Vlasov equation for /* and integrate by parts, this can be rewritten as 

c~ l A- 1 (d t V x B*)(x,t) 
= y j^(l + z®z) J (l--c- 2 v 2 )vd t f*{x + z,v,t)dv 

J^(l + z®z) J ((E* + c-'v x B*)f*)(x + z,v,t) dv 

+ 4c4 y y v2vd tf*( x + z ^' t ) dv - 

Therefore the solution E* of (|4.2|) has the representation 

E*(x,t) = 4vrA- 1 (Vp*)(x,t) + c- 1 A- 1 (^V x 
= -J\z\- 2 z P *(x + z,t)dz 

+^ J j \z\- 2 z{3{z-v) 2 -v 2 )(l-~c- 2 v 2 '}r{x + z,v,t)dvdz 
~2c2 / / l-^r 1 (1 + ^ «> ^) J {{E* + c~ l vx B*)f*){x + z,v,t) dvdz 

J J \z\~ 1 v 2 (v + (z -v)z)d t f*(x + z,v,t)dv dz. (4.3) 

Comparison with (jVPJ) and ()1.7|) reveals the analogy to E D at the relevant orders of c -1 . In 
particular, if we evaluate this relation at t — 0, the Banach fixed point theorem applied in C(,(1R 3 ) 
shows that for c > c* sufficiently large the function E*(x, 0) is uniquely determined by f°(x,v) = 
f*(x,v,0). Thus /° alone already fixes E° and 5°. Concerning the local and uniform (in c) 
existence of a solution to (jDVMc)) and the conservation of energy, one can use (|4.1|) and ()4.3|) to 
follow the usual method by setting up an iteration scheme for which convergence can be verified 
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on a small time interval; cf. Sect. 5.8]. Finally, by similar arguments as used in the proof of 
Theorem 1 1.1 1 it can be shown that solutions of (jDVMc|) approximate solutions of (jRVMc|) up to an 
error of order c~ 3 . □ 



5 Appendix 

5.1 Representation Formulas 

5.1.1 Representation of the approximation fields E D and B D 

Here we will derive the representation formula ()3.1|) for the approximate field E D from (jl.8|) . Since 
the calculations for the electric and the magnetic field are quite similar, we will only analyze in 
detail the electric field and simply state the result for its magnetic counterpart. 
From (jl.8)l we recall E D = Eq + c~ 2 E2, where 



E (x,t) = — J \z\ 2 z po(x + z,t) dz, (5.1) 

E 2 (x,t) = ~^J z d 2 po(x + z,t) dz — J Iz^dtjoix + z,t) dz — J \z\~ 2 z p 2 {x + z,t) dz, (5.2) 

cf. (IVPJ) and (|1.6|) . We split the domain of integration in {\z\ > ct} and {\z\ < ct}, and to handle 
the interior part {\z\ < ct} we expand the densities w.r.t. t about the retarded time 



t(z) :=t-c~ 1 \z\. 



To begin with, we have 



z\ z p (x + z,t) dz = — / \z\ z p (x + z,t(z)) dz 

z\<ct J \z\<ct 

-c _1 / M" 1 ^ d t po(x + z, i(z)) dz — c~ 2 - [ z d 2 po(x + z,i(z)) dz 



z\<ct ^ J \z\<ct 

-2- I u ^2n3, 



\z\ z (t — s) dtPo(x + z,s)dsdz. (5.3) 

z\<ct Jt(z) 



Using (|2.1|) and ()2.2|) . the last term is O cpt (c~ 3 ); note that \x\ < R for some R > together with 
the support properties of fo imply that we only have to integrate in z over a set which is uniformly 
bounded in c > 1. Since d t po + V ■ jo = by (jVP|) . we also find 

— c _1 / \A~ X z 9tPo( x + z i t( z )) dz = cT 1 / {z^zVx ■ jo(x + z,i(z)) dz 

J\z\<ct J\z\<ct 



-1 



c- 1 



z\<ct 



z\<ct 



\z\ x z v ■ V x fo(x + z, v, i(z)) dv dz 



\z\ 1 z v ■ (v z [fo(x + z, v, i(z))] + c 1 z d t fo(x + z, v, i(z)) j dv dz 



I + II, (5.4) 
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with 



c / \z\ z v ■ V z [fo(x + z,v,t(z))] dv dz 

J\z\<ct J 

J\z\<ct J ' ' 



+c (ct) / / z(z ■ v)f°(x + z,v) dv ds(z) 

J\z\=ctJ 

= — c -1 / / \z\~ 2 (v — 2(z ■ v)z)f (x + z,v,i(z)) dv dz 

J\z\<ct J 

+c _1 (ct)~ 1 / / z(z ■ v )f°(x + z, v ) dv ds(z); 

J\z\=ct J 

observe that i(z) = for \z\ = ct was used for the boundary term. Similarly, by (jVP 
II = c~ 2 / / \z\~ l z(z ■ v) d t fo(x + z, v, i(z)) dv dz 

J\z\<ct J 



c- 2 



/ \z\ 1 z(z ■ v)(v ■ VJo + E ■ V v f )(x + z,v,i(z)) dvdz 

J\z\<ct J 

V 2 • Ol -1 ^ ■ v)v ) i=lj2)3 f (x + z, v, i(z)) dv dz 



'\z\<ct ■ 

-c^ict)- 1 I I z(z-v) 2 f°(x + z,v)dvds(z) 



J\z\<ct J 

+c~ 2 / / • E f (x + z,v,i(z)) dvdz 

J\z\<ct J 

c~ 2 / / \z\~ 2 ((z ■ v)v + v 2 z — 3z(z ■ v) 2 ) f (x + z,v,i(z)) dv dz 

J\z\<ct J 



+c / \z\ z z ■ E p (x + z,t(z)) dz 

'\z\<ct 



-c" 2 (ct)- 1 / / z(z ■ v) 2 f°{x + z, v) dv ds{z) + O cpt {c~ 3 ). 

J\z\=ct J 

Next, due to (|2.1|) and (|2.2j) we also have 

]- i zd 2 p (x + z,i(z))dz = -c' 2 \ [ z d 2 p (x + z,t) dz + O cpt (c~ 

^ J\z\<ct * J\z\<ct 



-c- 2 - 
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Thus so far by (jSHJ) and Q - flHZZP , 

E (x,t) = — \z\~ 2 z p (x + z,t) dz — / |2|~ 2 2po(a ; + 2, dz 

J \z\>ct J\z\<ct 

\z\~ 2 z p (x + z, t) dz 



z\>ct 

\z\~ 2 z p (x + z,i(z)) dz — c~ 2 - J z d 2 p (x + z,t) dz 
— 2(z ■ v)z)f (x + z, v, t(z)) dv dz 

J\z\<ct J 

+ c- / / W -((, . v)v + A - • ,)») /„(* + „, £(,)) *, iz 

J\z\<ct J 
2 / l„|-l= = 



+c / |z| z z ■ E p (x + z,t(z)) dz 

\z\<ct 



+c (ct) I (z ■ v)z f°(x + z,v) dv ds(z) 

J\z\=ctJ 

-c^ict)- 1 [ f(z-v) 2 zf°(x + z,v)dvds(z) + O cpt (c- 3 ). 

J\z\=ct J 

Now we turn to E 2 , cf. (15. 2j) . Since i?2 enters with the factor c -2 , we first note that 

c~ 2 i / z d 2 p (x + z,t) dz 

^ J\z\<ct 

cancels a term on the right-hand side of (|5.8p . In addition, by analogous arguments, 
1^1 _ dtjo(x + z, t) dz 

/ |z|~ 1 9 t jo(a ; + z, i(z)) dz — c~ 2 / \ Z \~ X / d 2 jo(x + z,s)dsdz 

J\z\<ct J\z\<ct Ji(z) 



-c- 2 



z\<ct 

-c- 2 



'\z\<ct Jt(z) 

-c- 2 / / V 9tfo(x + z,v,t(z))dvdz + O cpt (c 3 ) 

\z\<ct 



/ \z\~ l / v(v V x f + E - V v f )(x + z,v,t(z))dvdz + O cpt (c- 3 ) 

J\z\<ct J 

M _1 vv ■ (v z [f (x + z,v,i(z))) + c^zdtfoix + z,v,i(z))j dvdz 
+cT 2 / l^l" 1 / vE -\7 v f (x + z,v,t(z))dvdz + O cpt (c- 3 ) 

J\z\<ct J 



c 2 / J \z\ 2 {z ■ v )v f (x + z, v, t{z)) dv dz 

J\z\<ct ■ 



+c 2 {ct) 1 / [z • v)vf°(x + z, v) dv ds(z) 

J\z\=ct J 

\z\- 1 (E oPo )(x + z,i(z)) dz + O cpt (c- 3 ). 



-c- 2 



\z\<ct 
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Finally, 

-c~ 2 / \z\~ 2 z p 2 (x + z,t) dz = -c~ 2 \z\~ 2 z p 2 (x + z,i(z)) dz + O cpt (c~ 3 ). (5.10) 

J\z\<ct J\z\<ct 



Therefore if we write 



c 2 E 2 (x,t) = c 2 - j z d 2 p${x + z,t) dz — c 2 f \z\ 1 d t jo(x + z, t) dz 

* J\z\>ct J\z\>ct 

—c~ 2 j \z\~ 2 z p 2 (x + z,t) dz + c~ 2 - j z d 2 p (x + z,t) dz 

J\z\>ct ^ J\z\<ct 

„-2 / IJ-lo • / | _ + \ j„ -2 / l_|-2; 



— c / \z\ d t jo(x + z,t) dz — c / 1^1 z p 2 (x + z,t) dz, 

J\z\<ct J\z\<ct 

use ()5.9j) and (|5.10p . and thereafter add the result to ([5.80 . it turns out that = + C 2 E 2 can 
be decomposed as claimed in ()3.1|) . 

Similar calculations for B D (x,t) = c~ 1 B 1 (x,t) using (jl.4|) yield 



B D (x,t) = c- 1 j 



\z\ 2 z x j (x + z, t) dz 



z\>ct 
1 / I I — 2 



+c / \z\ z x j (x + z,t(z)) dz 



\z\<ct 

— 2c~ 2 / / \z\~ 2 (z ■ v)(z x v) f (x + z,v,i(z)) dv dz 

J\z\<ct J 

+c~ 2 / M" 1 ^ x E p (x + z, i{z)) dz 

J\z\<ct 



-c- 2 (cty l {z-v){zxv)f°{x + z,v)dvds{z) + 0{c- 3 ). (5.11) 

J\z\=ct J 

5.1.2 Representation of the Maxwell fields E and B 

In this section we will verify the representation formula (|3.4j) for the full Maxwell field E, by 
expanding the respective expressions from (HJ |^] to higher orders. Once again the computation 
for the corresponding magnetic field B is very similar and therefore omitted. Let (/, E, B) be a 
C 1 -solution of ()RVMc|) with initial data (f°,E°,B°). We recall the following representation from 
ED (A13), (A14), (A3)]: 

E = Ed + Edt + Et + Eg, (5-12) 
B = B D + B DT + B T + B S , (5.13) 

where 

E D (x,t) = d t (— E° (x + ctu) duA + — i d t E(x + ctuj,0)duj, 
V 47r J\u\=i J 47r Jn=i 

E DT (x,t) = — (ct) _1 / / K DT (z,v)f°(x + z,v)dvds(z), 

J\z\=ctJ 

-2 



Ex(x,t) — — I \z\ / KT(z,v)f(x + z,v,t(z))dvdz, 

J\z\<ct J 

Eg(x,t) = ~c~ 2 \z\~ l Kg(z,v)(E + c^v x B)f(x + z,v,i(z)) dvdz, 

J\z\<ct J 
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and 



B D (x,t) = d t 



4tt 



B°(x + ctuj) duj 



|w|=l 



d t B(x + ciu;, 0) gL>, 



B DT (x,t) = (ci) 1 / / L DT (z,v)f°(x + z,v) dvds(z), 

J\z\=ct J 

B T {x,t) = cT 1 / \A~ 2 I L T (z,v)f(x + z,v,i(z))dvdz, 

J\z\<ct J 



B s (x 1 t) = c~ 2 l^l" 1 L s (z,v)(E + c^v x B)f(x + z,v,i(z)) dvdz, 

J\z\<ct J 

with z = \z\~ l z and t(z) = t — c~ l \z\. The kernels are given by 

K DT (z,v) = {l + c~ l z-vy l (z-c- 2 (z-v)v), 
K T (z,v) = (l + c- l z-v)- 2 {l-c- 2 v 2 )(z + c- l v), 

K s (z,v) = (l + c _1 £-£)~ 2 (l + c~V)~ 1/2 (l + c~ 1 z-v) + c~ 2 ((z-v)z-v)®v 



■(z + c 1 v) g) z 



i3x3 



and 



L DT (z,v) = (1 + c 1 z ■ v) 1 (z x c 1 v), 

L s (z, v) = (1 + c _1 z • v)' 2 (l + crV)~ 1/2 [(1 + c^z ■ v)z x (. . .) 

— c~ 2 (z x -0) <g) (cz + ■£>) 



i3x3 



Next we expand these fields in powers of c" 1 . According to ()2.5|) we can assume that the v-support 
of f(x,-,t) is uniformly bounded in x G M 3 and t G [0, T], say f(x,v,t) = for \v\ > P '■ = 
max ig [ 0i T] PvAi(t). Thus we may suppose that \v\ < P in each of the u-integrals, and hence also 
\v\ = (1 + c~ 2 v 2 )~ 1 / 2 \v I < \v\ < P uniformly in c. It follows that 

v= (l- l -c- 2 v 2 y + 0(c~ 4 ). 

For instance, for the kernel Kqt of E^t this yields 

K DT (z,v) = (l + c' 1 z-vy 1 (z-c~ 2 (z-v)v) 

= (l - c" 1 ^ • « + c" 2 (^ • w) 2 + C(c- 3 )) - cT 2 (^ • v)v + C(c" 4 )) 
= z - c~ l (z ■ v)z + tT 2 (z ■ v) 2 z - c~ 2 (z ■ v)v + 0(c~ 3 ). 

If we choose Ro > such that f°(x,v) = for |x| > i?o, then 



(ct)" 1 / / 0(c- i )l B{0 , Ro) (x + z)dvds(z) 

'\z\=ct J\v\<P 



Ct 



0(c 



1b(o,r ) + ctu) ds(u) \0(c ) 



w =1 
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by |2*T] Lemma 1], uniformly in x G M 3 , t G [0, T], and c > 1. Therefore we arrive at 
E DT (x,t) = -(ct)" 1 / / (z- c _1 (z-v)z + c' 2 [(z ■ vfz - (z-v)v])f°(x + z,v)dvds(z) 



\z\=ct 

,-3 



+0(cT d ). (5.14) 

Concerning Pr, we note that f(x, v, t) — for |x| > Po + PP =: Pi. Since, by distinguishing the 
cases \x — y\ > 1 and |x — y\ < 1, 

/ kT 2 lB(o,Hi)(a; + 2) cfe = / \z - y\~ 2 T-B^R^iy) dy = 0(1) 

J\z\<ct J\x-y\<Cb 

uniformly in x G M 3 , t G [0, T], and c > 1, similar computations as before show that 

E T (x,t) = - \z\~ 2 (z + c- 1 [v-2(z-v)z}+c- 2 [3(z-v) 2 z-v 2 z-2(z-v)v}) 

J\z\<ct J ^ ' 

f(x + z, v, i(z)) dv dz + 0(cT 3 ). (5.15) 
In the same manner, elementary calculations using also (j2.fi)) can be carried out to get 



E s {x,t) = -c 2 |^| 1 / [1 — z® z)(Ef)(x + z, v, i(z)) dvdz + C(c~ 3 ), (5.16) 

J\z\<ct J 

B DT (x,t) = (cty 1 I (c^z x v — c~ 2 (z ■ v)z x vj f°(x + z,v) dv ds(z) 

+0(c- 3 ), (5.17) 

B T (x,t) = cT 1 / \z\~ 2 hzxv-c~ 1 2vzzxv)f(x + z,v,i(z))dvdz + 0(c~ 3 ), (5.18) 

J\z\<ct J 



B s (x,t) = c~ 2 / \z\~ x / z x (Ef)(x + z,v,t(z))dvdz + 0(c- s ). (5.19) 

J\z\<ct J 

Next we consider the data term 

E D (x,t) = d t (— E°{x + ctu)duj J + — / d t E(x + ctu, 0) doj, 
V 471 " J\ u \=i / 4?r y|w|=i 

=: III + IV. (5.20) 
Since / 2 (z, v, 0) = by (JhVPj) . we have p 2 (x, 0) = 0. Thus we get from (jEJ), ffVPjt . and (fTHjl . 

P°(x) = P (x,0) + c~ 2 P 2 (x,0) 

= — y |2| _2 ^po(a; + 2, 0) cfc + c~ 2 J zd 2 p Q (x + z,Q) dz — J \z\~ l d t jo(x + z^) dz^j. 
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Using the formulas (|5.25j) . (|5.26|) . and (|5.24|) below, we calculate 

\z\- 2 z Po(x + ctu + z, 0) dz du 



|w|=l 



Pod/jO) / \y — x — ctu\ 3 (y — x — ctu) du dy 



= —An / \z\ 2 zpo(x + z, 0) dz, 

J\z\>ct 

/ z d 2 p (x + ctu + z, 0) dz du 

J\u\=lJ 

= d 2 p (y,0) / \y — x — ctu\~ 1 (y — x — ctu) du dy 

J J\u\=l 

= 4tt / (z - -(ct) 2 |z|~ 2 z) d 2 p (x + z, 0) dz + — ^ / 2 d 2 p (x + z, 0) dz, 
</|z|>ct v o / dcr ,/|«|<ci 

— / / l z l~ 1 ^io( a; + c ^ + z i 0) dz da; 

./U = l J 



dtjo{y,0) I \y-x-ctu\ l dudy 
M=i 



= — 47r / |z| 1 d t j (x + z, 0) dz / d t jo{x + z, 0) dz. 

J|z|>ci C ^ J\z\<ct 

Therefore we get 

III = d t (J^ J E°(x + ctu)du^j 

= d t ( —t \z\~ 2 z p (x + z, 0) dz + / zd 2 p (x + z, 0) dz 

V J\z\>ct ^ C J\z\>ct 

I \z\~ 2 zd 2 p (x + z, 0) dz + / zd 2 p (x + z, 0) dz 

" J\z\>ct J\ z \< ct 

— 4 [ \z\~ 1 dtjo(x + z, 0)dz - / % (^ + -2,0)dz J 
— — I \z\~ 2 z po(x + z, 0) dz H — c~ 2 / z d 2 p (x + z, 0) dz 

J\z\>ct 2 J\ Z \>ct 

— \t 2 I \z\~ 2 z d 2 p (x + z, 0) dz — c~ 2 / |z| _1 9 t jo(^ + z, 0) dz 

^ iui>ct J\z\>ct 



'|z|>ct J\z\>ct 

+ (ct)- 1 I zp (x + z,0)ds{z), (5.21) 

</ |z|=ct 

note that several terms have cancelled here. Now we discuss the second part IV of the data term 
E D , cf. ipT2Dj) . To begin with, by flQ, (fTIfl . and (jVFjl . 

S(x,0) = B°(x) = c _1 Bi(x,0) = c _i y |z|~ 2 z x j (x + z,0)dz 

= c~ x y y |z|~ 2 (z x v)fo(x + z, f , 0) dvdz, 
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due to jo — / vfo dv. Therefore using ()VP|) for f and integration by parts, we obtain 



Vx%0) = c _1 V x 




-i 






— c 




-i 




z\ 2 {z x v)f (x + z, v , 0) dv dz 
2 Vz/o x (z x v) dvdz 
~ 2 {(v ■ V x f )z - (z ■ V x f )vj dv dz 
- 2 z {-dtfo - E ■ V„/) dvdz 
z\~ 2 vz- V*[/ (. . .)] dv dz 
z\~ 2 z dtfo(x + z, v, 0) dv dz + c _1 47r / vfo(x,v,0)dv, 



by observing that in general for a continuous function g, 



\z\ 2 z ■ V ' z g(x + 2) dz 



Z\>£ 



g(x + eoS) duj — > 4iig(x), e — ► 0. 



w =1 



From Maxwell's equations we have dtE(x,t) = cV x B(x,t) — 4nj(x,t), and thus in view of 
j(x,0) = J vf(x, v, 0) dv = fvf°(x,v)dv = f vfo(x,v,0) dv, 

d t E(x, 0) = cV x 0) - 47rj(x, 0) 

|z|~ 2 z9 t / ( a; + -2; w ? 0) dt><iz + An / (n — v)f (x, v, 0) tfa;. 




Hence due to (15.251) , v — v = 0{c~ 2 ), and by [2U Lemma 1], 



IV 



4tt 



d t E(x + ctu, 0) du 



|oj| = 1 



47T 



|«|=1 




z\~ 2 z dtfo(x + ctu; + z, v, 0) civ duo 



+ c 1 (ct) / (v — v)f°(x + ctu,v) dv duj 

J\u\=l J 

-t I \z\' 2 zd t p (x + z,0)dz + O(c- 3 ). 

J\z\>ct 



(5.22) 



If we combine fl5~T2jl . (jOHjl . 4522), dS31, (fHUHJ) , and (|5T6|) . then we see that (jSJ) is 

satisfied. 

A similar calculation yields 

B D («,») = W W -«, **<* + *.<>)* + *"'*/' |,|-'**4ib(* + *,0)* 
</|z|>rf J\z\>ct 



' \z\>ct 

-c-\cty l 



z x Jo(^ + z, 0) g?s(z), 



(5.23) 



z\=ct 



and an analogous decomposition of B into £> = £> cxt + B int + B hd + O(o 
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5.2 Some explicit integrals 

We point out some formulas that have been used in the previous sections. For z G M 3 and r > 
an elementary calculation yields 

z - rcu = < , , , i - . (5.24) 

' 1 I 47rU : r < \z\ v 7 

Differentiation w.r.t. z gives 



\z — ruj\ 3 (z — ru) du = \ . , U ,_ 9 _ ' r > 2 _ (5.25) 
1 1 47r U \ z : r < \z\ v ' 



\u\=l 



Similarly, 

I z — ru I du 



'M=i 

and thus by differentiation 



Anr + ^-z 2 r 1 : r > \z\ 
4ir\z\ + ^r 2 ^! -1 : r < \z\ 



. z : r > \z\ 

3r 



z-ro;| \z-ru)du= I r . (5.26) 

'| w [=i I 47TZ — ^-r |z| z : r < |z| 

Finally, for zGl 3 \ {0} also 

\z - v\^\v \~ 3 vdv = 2t[z (5.27) 



can be computed. 

Acknowledgements: The authors are indebted to G. Rein, A. Rendall and H. Spohn for many 
discussions. 
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